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ABSTRACT: We study the correlation between the energy spectra of two disordered 
Hamiltonians of the form H a = Hq g + s a tp (a=l,2) with HQ a and (p drawn from random 
distributions. We calculate this correlation function explicitly and show that it has a 
simple universal form for a broad class of random distributions. 

a 

o 

In a series of recent papers [BZ1, BZ2, BZ3] (hereafter to be referred to as I, II, and III 
respectively) we studied the theory of random matrices [WIG,POR,MEH] and discovered 
a remarkably simple universality in the correlation between energy eigenvalues of random 
Hamiltonians. A brief review and summary of our work may be found in III. In particular, 
we find that for a large class of Hamiltonians, the smoothed connected correlation (to 
be defined below) between the density of energy eigenvalues at energies \i and v has the 
universal form 

smooth/ \ — 1 1 J (1 ~ X V) 



smootn ( x = - - J \ — \ . {universal) 

Pc ^ ' 2(iV 7 ra) 2 (x-y) 2 \ [(1 - x 2 )(l - j/ 2 )]V2 J 1 ' 

(where we have centered the spectrum at the origin.) Here a denotes the endpoint of 
the spectrum and we have introduced the scaled variables x = (i/a and y = v/a. We 
see that up to a factor of 1/a 2 , p^mooth j g e q U£ j ^ Q a universal function of x and y. The 
only dependence on the random distribution of Hamiltonians appears through a. This 



result is all the more remarkable since, by way of contrast, the density p(fj) is completely 
non-universal. This universality has also been obtained recently by Beenakker [bee] using 
another method. For p and v far from the endpoints of the spectrum, we see that in 
(universal) the factor in the curly bracket is essentially equal to one, and so the form 
the universal correlation function basically says that there is an algebraic long-distance 
correlation falling like l/(x — y) 2 . In III, we developed an efficient diagrammatic method 
to calculate the density of eigenvalues and the correlation between them. 

Meanwhile, in the recent literature Altschuler and collaborators [SSA] have considered 
the correlation between the energy spectra of two Hamiltonians H\ and H^. In particular, 
Simons and Altschuler considered the problem of a single electron moving in a ring threaded 
by a magnetic flux and with the electron scattering on impurities in the ring. The magnetic 
flux is then changed to some other value. Here H a = Hq q + ip for a = 1, 2, where HQ a are 
deterministic, representing the interaction of the electron with the magnetic flux, while tp 
is random, representing the interaction of the electron with the impurities. The correlation 
between the spectra of Hi and H2 is of interest in the physics of mesoscopic systems. In 
III, we determined this correlation for a certain class of distribution of <p. 

More generally, we may consider the correlation between the spectra of two Hamil- 
tonians H a = Hq g + s a <p for a = 1,2, where s a represent the strength of the disorder. 
Simons, Altschuler, and Lee [SSA] and more recently Narayan and Shastry [shastry] have 
considered this problem for the simpler case of Hqi = Hq2 = Hq. In this paper, we would 
like to demonstrate the power of the diagrammatic approach discussed in III by calculat- 
ing the correlation between the spectra of H\ and H2 for the more general case in which 
Hqi 7^ Hq2- Furthermore, instead of obtaining an implicit integral representation for the 
correlation, we will determine it explicitly. 

It may be useful to clarify the difference [diff] between our work and the work of 
Altschuler et al. They focus [SSA, shastry] on the correlation function p c (p, v) for p — v 
small, of order 1/N, (that is, for "distances" much smaller than the width of the spectrum), 
and indeed often for p and v both small. In III and in this work (and in the recent 
work of Beenakker [bee]) we look at the correlation function for p — v of order N®. In 
other words, we look at long distance correlation, while Altshuler and collaborators focus 
on short distance correlation. In fact, in I, by determining the large-index behavior of 
the relevant orthogonal polynomials, we were able to calculate explicitly the complete 
correlation function over all "distance" scales, from short distances through intermediate 
distances to long distances. The correlation function discussed in the literature appears as 
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a limiting case of the correlation function found in I. The work in I thus provides a bridge 
between the short distance discussion of Altschuler et al and the long distance discussion 
we gave in II and III. 

Let us begin by recalling some basic definitions and formulas. Consider an ensemble of 
N by N hermitean matrices <p defined by a probability distribution P(<p) and normalized 
by J dipP(ip) = 1. The large N limit will always be understood. Define as usual the Green's 
functions 

G <*»H>i^>i^) c (L4) 

where (0 (<£>)) = J d<p 0(<p)P(<p) and the subscript c indicates the connected Green's func- 
tion. The density of eigenvalues is then given by 

p(p) = {^tr6(n -<p)\ = ^ImG^ + ie) (1.6) 



and the correlation between eigenvalues, by 



(1.7) 

= (-1/4tt 2 )(G c (++) + G c (— ) - G c (+-) - G c (-+)) 

with the obvious notation G c (±, ±) = G c {ji ± ie, v ± iS) (signs uncorrelated). It is worth 
remarking that there are two limits involved: iV — > oo and e, 5 — > 0. As explained in III, 
we obtain the full correlation function (from which we can read off the correlation at all 
distances) by letting e, 5 — > first (the limit we took in I), and the smoothed correlation 
function by taking — > oo first. In this paper, we will look at the smoothed correlation. 

For applications to disordered system ip is to be thought of as the Hamiltonian. Its 
eigenvalues then describe the energy levels of the system. In some applications, (p is 
related to the transmission matrix [PICH]. From the density of and correlation between 
its eigenvalues we learn about the fluctuation of the conductance in disordered metals. 

In this work, we wish to study the correlation between the spectra of two different 
Hamiltonians Hi and H^- Evidently, we can determine this quantity by studying the 
generalized (connected) Green's function 

G(z,w) c = ( T7*r — - — Tffr — ) (crossG) 

\N z-tpiN w-(p2/ c 
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where to conform to field theoretic notation we have renamed H a as ip a , two "matrix scalar 
fields." Note that our starting point appears to correspond to one of the concluding results 
of Ref. shastry. (Cf. eq. (16) there.) 

For ease of presentation we will begin with the Gaussian distribution 

P(H 01 , #02, <p) = ie-TM^i+^^-HnV+^^oi-ffoa) 9 ]. (distribution) 

Zj 

Note that the factors of N are chosen in our definitions such that in the large N limit the 
density of eigenvalues has a finite (i.e., of order N®) width. For /j 2 — > oo, we recover the 
case Hqi = Hq2 considered earlier [SSA,shastry,BZ3]. In the opposite limit (i 2 = 0, the 
two Hamiltonians Hq\ and Hq 2 are completely independent. In evaluating G(z, w) c , we see 
that, since the integrand depends only on cp a , after using the definition cp a = Hq g + s a ip 
to eliminate Ho a and integrating over <p, we can effectively use P(<f) = ■^ 7 e^^ tr ^ aM ab L P h ^ 
with the 2-by-2 symmetric matrix M 2 defined by 

Mil = ( m l + z" 2 ) " M- 2 [(m\ + fi 2 )si - ^ 2 s 2 ] 2 , (ml) 

and similarly for Af| 2 , and 

M 12 = -V 2 ~ M- 2 [(m\ + v 2 )si - fj?s 2 ][(rnl + ^ 2 )s 2 - /Ai], (m2) 

with 

M 2 = m 2 + ^((m 2 + n 2 )s 2 a - 2siS2/U 2 - (m3) 

Let us now expand 

^ oo oo 

Gc(z, w) = j^J2H z m+l w n+l (tr^T(t)tr^(0)) c (expand) 
m=0 n=0 

Diagrammatically, we may borrow the terminology of large N QCD [thoo] and describe 
the expression for Gc(z, w) as two separate quark loops, of type z and type w respectively, 
interacting by emitting and absorbing gluons. A quark of type z can only emit and absorb 
gluons of "type 1" and a quark of type w and only emit and absorb gluons of "type 2." 
With a Gaussian distribution for cp a , a = 1,2, we can readily "Wick-contract" (expand). 
The inverse gluon propagator is given by 

(tPaijtPbkli®)) = j^ S il S jk^lb (glueprop) 

where the 2-by-2 matrix a 2 = M~ 2 . (For the sake of notational simplicity, we will define 
a 2 = (a 2 ) ii, t 2 = (cr 2 )22, and p 2 = (a 2 ) 12.) Thus, the gluon is represented by a double line 
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while a quark is represented by a single line. This convention greatly facilitates counting 
the powers of N. 

Let us begin by ignoring contractions within the same trace (in which case m and n 
are required to be equal). In the large N limit, the dominant graphs are given essentially 
by "ladder graphs" (with one crossing) which immediately sum to 

N 2 g c (z, w) = - — -2 -j— (ladder) 

V ZW' 



We next include Wick-contractions within the same trace in (tr^tr^) • Graphically 
these contractions correspond to decorating the ladder graphs by vertex and self energy 
corrections. Summing the vertex corrections, we see that we have to multiply the expression 
in (ladder) by two factors, the factor (1 — fr) 2 an d a similar factor with z and a replaced 
by w and r respectively. Finally, summing the self energy corrections we see that the bare 
quark propagator 1/z gets dressed to 

G(z, a) = (z — \/ z 2 — 4a 2 ) (green) 

(and similarly for l/w of course with a replaced by r.) Note that the bare propagators 
are recovered in the limit cr, r — > 0. Putting these various factors together, we obtain the 
remarkably compact result 

- 2 r G 2 (z,a) 1 [ G 2 (w,t) 



N 2 G c (z,w) P 



(l-p 2 G(z,a)G(w,r)) 2 



l-a 2 G 2 (z,a) 



1-t 2 G 2 (w,t) 

(connected2point) 

As indicated in III, we find it convenient to introduce angular variables: sin# = ji/2a and 
sin(/> = u/2t. As 9 and (p vary from —n/2 to tv/2, \i and v vary over of the width of their 
respective spectra. These angular variables are the natural scaling variables to use in this 
class of problems dealing with eigenvalues of random matrices. From (green), we have 

G(p + ie, a) = -ii]e^ 9 /a (cutGl) 

where r\ = the sign of e and 

G(v + id, t) = -iie l ^/r (cutG2) 

with £ = the sign of 5. A straightforward computation of the appropriate generalization 
of (1.7) then yields 

- IQti 2 N 2 Pc (^ v) 

1 J 1 + ch u cos(6> + 4>) 1 - ch u cos(6> - <j>) \ (correlation) 

ar cos 9 cos <fi \ [chu + cos(# + (p)] 2 [ch u — cos(# — (p)] 2 
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Here we have defined 



1 



2 \2i 



u = \og(ar/p') = - log[M^M| 2 /(Mf 2 ) 



(defineu) 



a rather involved function of the parameters m?, to 2 , to 2 , and p 2 that appear in the 
distribution (distribution). Perhaps remarkably, aside from an overall factor of or the 
correlation function, when expressed in terms of the appropriately scaled variables 9 and 
4> depends on the probability distribution only through this particular combination. This 
represents already a form of universality or scaling. 

The reader familiar with III would recognize that this is the same expression we ob- 
tained there (see eq. (2.18) in Ref. BZ3) for an apparently totally different problem, with a 
totally different definition of u. In III we studied the time-dependent correlation between 
the eigenvalues of time-dependent matrices <p(t) taken from a probability distribution de- 
fined by 



p (<fi) = 7F77 ex P 



dt Tr 



-T 



<pK(±)<p + V(<p) 



(timedep) 



with K any reasonable function and with T — > oo. Here "time" may correspond to some 
external parameter we are allowed to vary. For this problem, we defined e~ M W to be 
j ^.^t^ij x(oo)). For the problem considered in this paper, in contrast, e~ u is defined 
simply in terms of the widths of the distributions from which various random Hamiltonians 
are drawn. Evidently, we have in this way found a "mapping" between two apparently 
unrelated problems involving disordered Hamiltonians. A particularly simple case occurs 
when K(-i^) = —-^z in which case u(t) is essentially time. It is in this sense that loosely, 
one may speak of u as defined in (defineu) as "time," even though there is no notion of 
time as such in the problem studied in this paper. 

Various limits can now be studied. In the large "time" limit u — > oo, we obtain 



4n 2 N 2 p c (p, is) — > e u tan#tan ( 



(longtime) 



For 6 = <fi (note that this does not mean ji = v when a is not equal to r) and small "time" 
u we find 



8n 2 N 2 p c (p, v) = 



1 1 

or cos 2 9 u 2 



+ 



For 9 7^ 4> an d u « 9 — 4> we obtain 



16tt 2 N 2 Pc (p, v) = 



or cos 9 cos 4> I 1 + cos(6' + (f>) 



1-2 cos(6> + 4>) 2 
1 + cos(6> + 4>) U 



(samepoint) 

+ ((j) -> -(f) + ?r) 
(smalltime) 



We can also study various regions of parameter space, of course. Consider (A) \i — > oo 
so that Hqi = Hq2 = Hq in which case it is convenient to define m 2 , = m 2 + m|. Case 
(Al): with mo — > oo (so that i^o = 0) and s\ = S2 we recover the universality (universal) 
found in I and in Ref. bee. Case (A2): with rriQ » m so that the disorder is large, 
we find that e u = ar / p 2 = 1 and thus the same universal correlation as in case (Al) 
suitably scaled by an overall factor of err, but here there is no requirement on s a - Case 
(A3): with S2 equal to zero, so that we are studying the correlation between the spectra 
of a "bare" Hamiltonian Hq and a Hamiltonian disturbed by an external perturbation 
Hq + si(p, we find that e u = 1 + m^s^/rn? . Note that no general statement on how "time" 
t maps onto the strength si of the external perturbation can be made, but in the simplest 
case of identifying "time" mentioned above we have t oc logsi for large strength of the 
perturbation. 

Another interesting class (B) is defined by the limit \i — > so that Hq\ and Hq2 

c\ O c\ c\ 

are not correlated at all. After some drastic simplification, we find a = + sfm~ , 
t = m 2 + s^mT , and p = siS2fn~ . As a check on the formalism, we can consider 
some special cases. For instance, with m a — > oo, a = 1, 2, Hqi and Hq2 are set to zero, 
in which case we recover essentially our previous result (universal). As another example, 
consider the analog of case (A3): with S2 = we obtain u — > oo and as expected p c = 
since now Hi and H2 are no longer "linked" by ip. 

All these results are derived with the Gaussian distribution (distribution). We now 
remark on how these results may hold for more general distribution. In III we identified 
two classes of random matrix distribution which we refer to as the Wigner class and 
the trace class. In the Wigner class, we consider an ensemble of matrices whose matrix 
elements are independently distributed according to some probability distribution (the 
same distribution for all the elements.) Let us focus on the example in which the probability 
of the distribution matrix element cpij is given by 

P(ip tj ) oc e _Ar2(l ^' |2_ ^ )2 . (wigner) 

(Here ip represents any of the matrices in (distribution) generically.) It is easy to see that 
in the Feynman diagrams involving the quartic interaction ~ A^ 2 |</7^| 4 the indices % and 
j are tied together so to speak and thus these graphs are suppressed by a power of A^ 1 
relative to the graphs in which this interaction vertex does not appear. Reasoning along 
this line, we see immediately that the results of this paper have a generality far beyond 
the Gaussian distribution. Our results are universal. This universality is much harder to 
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prove for the so-called trace class in which the generic random matrix ip is taken from a 
distribution of the form P((p) = ^ e ^ NtrV ^ . We conjecture however that our results may 
also hold for random Hamiltonians defined with this class of probability distribution. In 
II, we had shown that the universal correlation in (universal) indeed holds for this trace 
class and for a much broader class generalizing the trace class. In III, we also calculated 
the correlation when Hq\ and Hq2 are set to equal deterministic, rather than random, 
Hamiltonians W\ and W% say. This case may be considered in the present context by 
adding to the logarithm of the distribution in (distribution) terms like —a^(HQ a — W a ) 2 
and taking a a — > oo. The gluons ip\ and <p2 are now coupled also to external sources. 

In conclusion, we have extended the universality in (universal) to a much broader 
class of problems involving random Hamiltonians. Expressed in terms of suitably scaled 
variables, the correlation between energy eigenvalues is universal, certainly in the Wigner 
class, and most likely in the trace class as well. The case Hq\ = Hq2 discussed previously 
in the literature can be recovered as a special case within our discussion. 
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